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Summary 

Design sensitivity is central to most optimization methods. The analytical sensitivity expression for 
an indeterminate structural design optimization problem can be factored into a simple determinate term 
and a complicated indeterminate component. Sensitivity can be approximated by retaining only the 
determinate term and setting the indeterminate factor to zero. The optimum solution is reached with the 
approximate sensitivity. The central processing unit (CPU) time to solution is substantially reduced. The 
benefit that accrues from using the approximate sensitivity is quantified by solving a set of problems in a 
controlled environment. Each problem is solved twice: first using the closed-form sensitivity expression, 
then using the approximation. The problem solutions use the CometBoards testbed as the optimization 
tool with the integrated force method as the analyzer. The modification that may be required, to use the 
stiffener method as the analysis tool in optimization, is discussed. The design optimization problem of an 
indeterminate structure contains many dependent constraints because of the implicit relationship between 
stresses, as well as the relationship between the stresses and displacements. The design optimization 
process can become problematic because the implicit relationship reduces the rank of the sensitivity 
matrix. The proposed approximation restores the full ra nk and enhances the robustness of the design 
optimization method. 


Introduction 

Design sensitivity is central to most optimization methods. It can be a major contributor to the 
number of calculations in optimization. The computation of efficient design sensitivity for structural 
problems has drawn considerable attention. NASA organized a conference on the subject matter (ref. 1). 
It is also well documented in the literature (refs. 2 to 6). General-purpose codes provide for the 
calculation of sensitivity for stress and displacement constraints (ref. 7). The automatic differentiation 
of the FORTRAN code, ADIFOR (ref. 8), has also been suggested to calculate sensitivity. In such a 
circumstance, we ask and attempt to answer a question about the precision of the sensitivity of stress and 
displacement constraints in the design optimization of an indeterminate structure: “Is the optimization 
process robust when the design sensitivity matrix is highly accurate?” The contrary may be true. The 
performance of an optimization method can be improved when the analytical sensitivity is replaced by a 
determinate approximation. The approximate sensitivity matrix is not only adequate, but it should be 
preferred in the design calculation of an indeterminate structure. Optimization, in other words, requires 
sensitivity, but approximate gradients are quite satisfactory. To illustrate the benefits that accrue from the 
approximations, the authors used several indeterminate trusses as numerical examples because design 
optimizations have been completed for such structures. The concept, however, should be extendable to 
other types of structures, such as beams, framework, and shell structures. 

Consider a truss that is made of n bars with r dependent members. The n bar areas are treated as the 
design variables. Approximate sensitivity works well because of three attributes special to an 
indeterminate truss. 
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(1) Dependent stresses: In an indeterminate truss, r out of n bar stresses {a} are dependent. Stresses 
are dependent because of the r compatibility conditions, which can be written as [c] {a} = {0} . The 

r x n sparse matrix is independent of the design variables. 

(2) Dependency of stresses and displacements: The m = n - r number of displacements {X} are 
dependent on the n bar stresses: {X} = [i?] {a}. The m x n sparse matrix is independent of the 
design variables. 

(3) Active constraints: In structural design, the number of active constraints can exceed the number 
of design variables. In an optimization algorithm, the singularity condition can be alleviated by restricting 
the number of active constraints to not exceed the number of design variables. 

In an optimization algorithm, the calculation of the search direction {d} requires the constraint 
gradient matrix [Vg]. An example of the use of the sensitivity matrix to generate the direction follows: 


M=[e]{v/} 

a) 

r n M-[vg][[Vgf[vg]l _1 [vgf 

[ 2 ] = 1 


Here /is the objective function, and 


[tfH'Hvsfef M Vsf 


The matrix 


[Vgf [Vg] 


becomes singular for each of the three special attributes, items (1) to (3), 


stated above. An optimization algorithm may yield a solution despite the singularity condition because the 
[Q\ matrix is approximated most often; it is seldom calculated in closed form. It is reinitialized into an 
identity matrix when corruption is suspected. The proposed approximate sensitivity of stress and 
displacement constraints alleviates the singularity condition in the design optimization of an indeter- 
minate truss. The solution is reached with fewer calculations because the optimization process becomes 
more robust, and the sensitivity is generated with a trivial amount of computations. The benefit that 
accrues from the use of approximate sensitivity is shown through the solution of a set of problems that 
were selected from the literature. Each problem was solved twice in a controlled environment, first using 
the closed-form gradient, then with an approximation. A comparison of the two optimum solutions 
quantified the benefit. The underlying cause of the benefit was investigated through a discussion of the 
nature of structural design optimization problems. 

This paper is organized into six subsequent sections. The design problem is formulated in the second 
section. The analysis and optimization tools are discussed in the third and fourth sections, respectively. 
Solutions to a set of problems are given in the fifth section, followed by a discussion and conclusions in 
the sixth section. A symbols list is given in the appendix to aid the reader. 


Design Optimization Problem 

Minimum weight is the objective of the truss design problem. The bar areas A, are considered as the 
design variables. Limitations specified on the bar stress a, and nodal displacement uj form the behavior 
constraints. The design optimization is cast as the following mathematical programming problem. 
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Find n design variables, here bar areas j A j : 

n 

To minimize the weight of the truss, W (A) -I p iMi 

7=1 


under n stress constraints, 


gM) = 




ho 


- 1<0 


( 2 ) 


and 777 displacement constraints, 


g n +j( A ) = 


u j 


u j o 


- 1<0 


Here p, and £ t are the weight density and length, a, and a, 0 are the stress and allowable strength, and Uj 

and 77/o are the nodal displacement and limitation, respectively. For a large problem, the number of design 
variables can be reduced by linking bar areas. Likewise, a small number of critical constraints can be 
separated and used in the design calculations (ref. 9). 


Design Update Formula 

A key formula to update the design variables (here area A,) in a nonlinear programming algorithm at a 
k' h intermediate iteration can be written as 


{ A )k { A )k-\ + a £-l WIat-1 


(3) 


The step length a*-i is calculated to minimize the weight along the direction {d} k _\ inside the feasible 
domain. All n + m constraints should be used to define the feasible space. The sensitivities of the set of 
active stress and displacement constraints are required to calculate the direction vector The quality 
of the direction vector is dependent on the accuracy of the sensitivity matrix. A spurious direction would 
be generated if the sensitivity matrix was rank deficient. 

Consider the j th stress constraint. Its closed- form gradient can be expressed as the sum of two factors: 



determinate 


+ 



indeterminate 


(4) 


The gradient expression given by equation (4) has to be adjusted for the absolute value in the 
constraint, which however, poses no limitation to the discussion here. The first factor { Va ; } dcle """ iate 
is applicable to determinate as well as indeterminate trusses. This vector has only one nonzero entry, 
which is the negative ratio of the member force to the square of the bar area ( -FI A 2 ). The second term 
{Vo / } mdel ™ mat£ accounts for the effect of indeterminacy. It is not a negligible factor. It can be fully 
populated, and its calculation is computationally intensive. The proposition is to drop the second term 
{Va / }' ndctcrminate in design optimization even when it is nontrivial. The nature of the gradient of the 
displacement constraint is quite similar to that of the stress constraint. Again, the proposition is to retain 
only the simple determinate factor. 


Analysis Tool 

An analysis tool is required to calculate the stress and displacement constraints and their sensitivities. 
Here, the integrated force method (IFM) (ref. 10) is employed. The structure of the IFM equation is 
suitable to calculate the closed-form sensitivities because the sizing design variables of a structure (here 
bar areas) are retained in a pristine state in the concatenated flexibility matrix [G], In addition, IFM has 
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two distinct sets of equations. The first set, with internal force as the primary unknown, is differentiated 
to obtain the sensitivity of stress. Likewise, the sensitivity of displacement is recovered by differentiating 
the second set of equations. The adjustment that may be required for the stiffness method of analysis is 
also discussed in this paper. The 1FM equations to calculate forces and back-calculate displacements are 
as follows: 

Internal forces {F} are calculated from the governing 1FM equation: 


[S]{F}={P} 


Displacements {X} are back-calculated from the forces: 


Flere 


{X}=p][G]{F} 



[C][C] 


and 


[5] 

[B] 

[Q 

[G] 


{P} 


in 

{ 5 R} 

w 

[J] 


n x n governing matrix 
m x n equilibrium matrix 
r x n compatibility matrix 
n x 77 flexibility matrix 
77-component force vector 

r p 


77-component load vector, {P} = < 


> 


S R 


777-component mechanical load vector 
r-component initial load vector, r=n — m 
777-component displacement vector 
first 777 rows of [[5] '^matrix of dimension m x n 


(5) 

(6) 


The sensitivity matrix for the stress and displacement constraints for an 77 -bar truss with r dependent 
members is obtained by differentiating the IFM equations. The closed-form sensitivity matrix for stress 
has the following form (ref. 11): 



- 

r • -T 

- detenninate 

- 

r- -i-l 

[Vg] = [Vgj, Va 2 ,..., Va„ ] = 

- 

• F 

+ 

[d] t 

' 1 

A-. 


-1 indeterminate 


(7) 


The expression given by equation (7) should be adjusted for the allowable strength prior to its use in 
design optimization. The rank of the ?7 x 77 sensitivity matrix [Vg] in equation (7) is reduced to m = n - r 
when both terms are retained. The recommendation is to use only the first term in equation (7), which is 
superscripted “determinate.” It is a diagonal matrix with full ra nk n. The proposition is to drop the second 
term that is superscripted “indeterminate.” The calculation of this term is computation intensive. 
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The closed-form displacement sensitivity follows: 


[vx] = 




determinate 


+ 




indeterminate 


( 8 ) 


The proposition is to use the first factor with superscript “determinate” in design optimization. The 
definitions of the symbols in equations (7) and (8) are as follows: 


l D H s ] 


-i 


[0] 


[C] 


[/] 


The elements of the diagonal matrix ' -Sjg - . J arc given by 

SiiFi liFi 




4 4 E i 


• J_ 
A- 


F 

~A*-- 


G . 


, and 'F ■. J are diagonal matrices of dimension n x n whose elements 


IF t- 

are — , — — 4 — , and , respectively. 

Ai Af AfEi 

The determinate factor in the displacement sensitivity can be specialized for an n-bar indeterminate 
truss as 


[VX] 


detenninate 


mM 1 


0 


-FI 

a 2 e 




(9) 


The displacement sensitivity given by equation (9) is a function of the bar length t. Young’s modulus 

E, and areas A because displacement is a global variable. Calculation of the determinate sensitivity for the 
displacement essentially requires a back-substitution step with the factored form of the inverse of the S 
matrix. It is important to observe the similarities and differences in the sensitivity expressions of the stress 
and displacement. 

(1) Sensitivities of both the stress and displacement contain the member forces {F\ and the square of 
bar areas {A Even for an active displacement constraint, the design in essence is modified through the 
member force. 

(2) The geometry of the truss is not explicitly contained in the sensitivity expression for stress 
because it is a local variable. 

(3) Because displacement is a global variable, its sensitivity expression explicitly contains the 
geometrical or configuration parameters, the material property, and the design variables. 
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Optimization Tool 

The consequences of using approximate sensitivity in design optimization are demonstrated through 
the solution of a set of problems. Solutions were generated within the framework of the design 
optimization testbed CometBoards (ref. 12). Each problem was solved twice, first using the determinate 
sensitivity, then with the full closed-form expression. Problems were solved in a controlled environment 
on an SGI workstation running the IRIX 6.5 operating system (Silicon Graphics, Inc., Mountain View, 
CA). Identical convergence and stop criteria were used for the optimization algorithm. For large 
problems, we reduced the number of design variables and constraints by utilizing the design variable 
linking and constraint formulation features available in CometBoards. A sequential quadratic program- 
ming algorithm, SQP (ref. 12), was the primary optimizer. This algorithm was supplemented, when 
required, by a modified method of feasible directions (mFD) and a sequential unconstrained minimization 
technique (SUMT). 

Research to compare different optimization algorithms and alternate analysis methods for structural 
design applications has grown into a multidisciplinary design testbed that is still referred to by its original 
acronym, CometBoards, which stands for COMparative Evaluation Test Bed of Optimization and 
Analysis Routines for the Design of Structures. The modular organization of CometBoards, shown in 
figure 1, allows for quick testing of innovative methods (or computer codes) in a controlled environment 
through its soft-coupling feature. Optimizers and analyzers are two important modules of CometBoards. 
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The optimizer module includes a number of algorithms: the fully utilized design method, optimality 
criteria methods, the method of feasible directions, mFD, three different versions of SQP techniques, 
SUMT, the sequence of linear programming, a reduced gradient method, and others. Likewise, the 
analyzer module includes several structural analysis codes, an aircraft flight optimization analyzer, a jet 
engine performance program, and others. CometBoards has several unique features, including a multiple 
optimizer cascade strategy, design variable and constraint formulations, a global scaling strategy, analysis 
and sensitivity approximations, regression and neural network approximators, and substructure 
optimization in sequential as well as in parallel computational platforms. CometBoards can accommodate 
up to 10 different disciplines, each of which can have a maximum of five subproblems. The testbed can 
optimize a large system that can be defined in as many as 50 different subproblems. Alternatively, a 
component of a large system can be optimized. 

Numerical Examples 

Solutions were generated for a set of six examples. The number of design variables ranged between 
3 and 23 linked variables. The constraints ranged between 7 and 312. First, we summarize each problem 
and provide its optimum solution. This is followed by the CPU time to solution. We conclude the section 
with a discussion on the search direction. 

Numerical Example 1: A Three-Bar Truss 

The optimum solution was calculated for the three-bar steel truss shown in figure 2 using the 
determinate term as well as the full-sensitivity expression. The truss in this figure was subjected to three 
different load cases. A 100-kip load along the negative incoordinate direction at the free node 1 was the 
first load case. The second and third load cases consisted of a 100-kip load in the positive and negative 
x-coordinate directions, respectively. The three bar areas were the design variables, and minimum weight 
was the objective. The allowable stress was 20 ksi for each member. The displacement limitations were 
0.25 and 0.50 in. at node 1 along the x and y directions, respectively. There were a total of nine stress and 
six displacement constraints. The problem was solved using the SQP algorithm, and the optimum solution 
is given in table 1. 



Load, 

kip 

Load condition 

1 

2 

3 

Px 

0 

100 

0 

Py 

-100 

0 

-100 
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TABLE I.— OPTIMUM SOLUTION FOR THE THREE-BAR TRUSS 


Weight, 

lbf 

Design variables, 
in. 2 

Active constraints 

Sensitivity 


^1 

a 2 

^3 

Stress 

Displacement 


163.6 

1.89 

0.33 

1.89 

1 

2 

Detenninate 

163.8 

1.89 

0.33 

1.89 

1 

2 

Analytical 


The SQP algorithm converged to the same optimum solution for the determinate as well as for the 
full-sensitivity expression. At the optimum, the rank of the sensitivity matrix was 2 and 3 for the 
analytical and determinate sensitivity expressions, respectively. The three-variable problem had three 
active constraints. 


Numerical Example 2: A Tapered Five-Bar Truss 

The optimum solution was obtained for the tapered five-bar steel truss shown in figure 3. The truss 
geometry and loads are depicted in the figure. The allowable stress was 20 ksi, and the displacement 
limitations were 0.25 and 0.50 in. at node 2 along thex and y directions, respectively. The five-bar areas 
were the design variables, with a 0.25-in. 2 lower bound on bar areas. The problem had a total of seven 
constraints. It was solved using three algorithms: SQP, mFD, and SUMT. The optimum solutions are 
given in table II. 


y 
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TABLE II.— OPTIMUM SOLUTIONS FOR THE FIVE-BAR TRUSS 


Method 

Weight, 

lbf 

Design variables, 
in. 2 

Active constraints 

Sensitivity 

^1 

^2 

^3 

A 4 

^5 

Stress 

Displacement 

SQP 

645.0 

2.21 

1.41 

1.65 

1.60 

0.25 

5 

2 

Determinate 

mFD 

646.7 

0.84 

2.78 

0.25 

2.98 

0.63 

4 

2 

Determinate 

SUMT 

646.8 

1.74 

1.89 

1.17 

2.09 

0.34 

5 

2 

Determinate 

SQP 

646.2 

2.48 

1.15 

1.93 

1.34 

0.25 

4 

2 

Analytical 

mFD 

648.5 

1.81 

1.83 

1.24 

2.03 

0.37 

4 

2 

Analytical 

SUMT 

647.7 

1.81 

1.82 

1.25 

2.01 

0.36 

4 

2 

Analytical 


This example appears to have multiple optimum solutions with about the same weight of 646 lbf. The 
SQP, mFD, and SUMT algorithms converged to different optimum solutions with a small variation in the 
weight. For determinate sensitivity, the mean weight was 646.2 lbf, with a maximum deviation of about 
0.2 percent. For the closed-form sensitivity, the weight was marginally higher at 647.5 lbf, also with a 
0.2-percent variation. The variation in area was rather wide. Consider, for example, the area for the first 
bar. For the full-sensitivity expression, its mean was 2.0 in. 2 , with a maximum deviation of 22 percent. 
For the approximation, the mean was 1.6 in. 2 , with a maximum variation of 47 percent. The sum of the 
areas for bar 1 and bar 2 remained at 3.62 in. 2 When bar 1 became heavier, bar 2 became lighter and vice 
versa. The active constraint set included two displacement limitations and either four or five stress 
constraints. The determinate sensitivity performed at about the same level as the closed- form gradients. 
The five-variable problem with six active constraints was not a well-posed mathematical programming 
problem (ref. 9). 


Numerical Example 3: A Forward-Swept Wing 

A forward-swept wing made of aluminum was modeled as a space truss with 135 bars, as shown in 
figure 4. It was subjected to loads at the wing tip that induced flexure in the x-z plane and torsion in the 
y-z plane. The allowable stress was ao = 10 ksi for all members. Displacements along the z-coordinate 
direction were constrained at nodes 10 and 30 with a 2-in. limitation. The 135 bar areas were grouped to 
obtain 23 linked design variables. The problem had 135 stress and two displacement constraints. The 
optimum solutions obtained by the SQP and SUMT algorithms are given in table 111. The SQP algorithm 
converged to the minimum weights of 4218.5 and 4217.9 lbf for analytical and determinate sensitivities, 
respectively. For the SUMT algorithm, the weights were 4072.7 and 4071.5 lbf, respectively. The 
3.5-percent difference in the minimum weight between the algorithms could be attributed to the 
complexity of the problem. The mean values of the design variables are depicted in table 111. Eight 
variables converged to the lower bound (0.25 in. 2 ). There were six and nine variables, above and below 
the mean value, respectively. A consistent set of active constraints was generated for the determinate as 
well as the analytical sensitivities. However, there were more active constraints than the number of design 
variables. The mFD algorithm converged to a heavy design that was considered incorrect and excluded 
from discussion. The solution with approximate sensitivity required a much smaller number of 
calculations. The central processing unit (CPU) time to solution for the determinate and the closed-form 
sensitivities was reduced by factors of 7.8 and 5.5 for the SQP and SUMT methods, respectively. 
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TABLE III.— OPTIMUM SOLUTIONS FOR THE FORWARD-SWEPT WING 


Method 

Weight, 

lbf 

Mean value of 
design 
variables, 
in. 2 

Active constraints 

CPU 

time, 

sec 

Sensitivity 

Stress 

Displacement 

SQP 

4217.9 

7.34 

20 

1 

8.7 

Determinate 

SUMT 

4071.5 

7.11 

24 

1 

6.2 

Determinate 

SQP 

4218.5 

7.35 

20 

1 

67.5 

Analytical 

SUMT 

4072.7 

7.11 

24 

1 

34.0 

Analytical 


Numerical Example 4: A Trussed Ring 

The design of the trussed steel ring shown in figure 5 was considered next. The ring was made of 
60 bars and had inner and outer diameters of 180 and 200 in., respectively. It was fully restrained at 
node 10 and free to move only along they direction at the diametrically opposite node 16. The ring was 
subjected to two load conditions. The first load condition consisted of a 40-kip compression along the 
ring’s horizontal diameter, which was applied at nodes 1 and 7. In the second case, a 40-kip load was 
applied at node 4 to induce compression along the vertical diameter. 
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40 kip 

(load condition 2) 


m 



Figure 5. — Sixty-bar trussed ring. 


The 60 bar areas of the truss were grouped to obtain 16 linked design variables. The ring had 
60 stress constraints (with a yield strength of 20 ksi) for each load condition. The distortions of the 
ring along the horizontal and vertical diameters were controlled through a 4-in. displacement limita- 
tion specified at nodes 1, 4, and 7 for each load condition. The problem had a total of 120 stress and 
6 displacement constraints. Optimum solutions generated by SQP and SUMT algorithms are given in 
table IV. Mean, maximum, and minimum values are given for the 16 design variables. 


TABLE IV.— OPTIMUM SOLUTION FOR THE TRUSSED RING 


Method 

Weight, 

lbf 

Design variables, 
in. 2 

Active constraints 

CPU 

time, 

sec 

Sensitivity 

Mean 

value 

Variation 

Stress 

Displacement 

Minimum 

Maximum 

SQP 

799.9 

3.15 

2.18 

4.05 

28 

1 

2.3 

Determinate 

SUMT 

797.7 

3.15 

2.18 

4.01 

28 

1 

1.3 

Determinate 

SQP 

799.9 

3.15 

2.18 

4.04 

28 

1 

7.4 

Analytical 

SUMT 

798.0 

3.15 

2.19 

4.00 

28 

1 

7.1 

Analytical 


Both SQP and SUMT algorithms with approximate as well as full closed-form sensitivities converged 
to the same solution with a 0.25-percent deviation in the minimum weight. The CPU time to solution was 
321 and 318 percent faster for approximate sensitivity with the SQP and SUMT algorithms, respectively. 
The 16-variable problem had 29 active constraints. 
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Numerical Example 5: A 25-Bar Truss 

The aluminum tower modeled as the 25 -bar space truss shown in figure 6 was designed for minimum 
weight under stress and displacement constraints. It was subjected to two load conditions. In the first load 
case, node 1 was subjected to three load components: 5-, 20-, and -5 -kip forces along the x, y, and z axes, 
respectively. In the second load case, node 2 was subjected to -20- and -5-kip forces along the y and 
z axes, respectively. 

The allowable stress was 1 0 ksi, and the displacement limitation was 1 in. in all three directions 
for the six free nodes. The 25 bar areas were linked to obtain eight design variables. The problem had 
25 stress and 18 displacement constraints. The optimum solutions calculated by the SQP and mF D 
algorithms are given in table V. The mean value and variance are given for the eight design variables. The 
eight-variable problem had nine active constraints. Both the SQP and mFD algorithms converged to the 
same optimum solution for the determinate as well as the analytical gradients. 
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TABLE V.— OPTIMUM SOLUTIONS FOR THE 25-BAR TRUSS 


Method 

Weight, 

lbf 

Design variables, 
in. 2 

Active constraints 

Sensitivity 

Mean value 

Variance 

Stress 

Displacement 

SQP 

190.5 

0.49 

0.13 

7 

2 

Determinate 

mFD 

190.6 

0.49 

0.13 

7 

2 

Determinate 

SQP 

190.5 

0.49 

0.13 

7 

2 

Analytical 

mFD 

190.4 

0.49 

0.13 

7 

2 

Analytical 


Numerical Example 6: A 20-Bay Truss 

The minimum-weight design was calculated for the 20-bay steel truss shown in figure 7. The 
structure in the figure was subjected to three load cases. The first load case consisted of forces in the 
negative y-coordinate direction along the bottom chord nodes: -40 kip at the midspan and -1 kip at the 
other nodes. For the second load case, all the top chord nodes were subjected to a 3-kip force along the 
x-coordinate direction. For the third load case, all the bottom chord nodes were subjected to a -3-kip force 
along the negative x-coordinate direction. 

The allowable stress was a 0 = 20 ksi. A displacement limitation of 0.5 in. was imposed at the midspan 
nodes 2 1 and 22 along the x- and y-coordinate directions, respectively. The 101 bar areas were grouped to 
obtain five linked design variables. The first two design variables represented the bar area of the top and 
bottom chord members, respectively. All 21 vertical bar areas were grouped to obtain the third variable. 
The last two design variables represented the bar area of the leading and lagging diagonal members, 
respectively. The five-variable problem had a total of 312 stress and displacement constraints. Optimum 
solutions obtained by the SQP algorithm are given in table VI. 
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TABLE VI.— OPTIMUM SOLUTIONS FOR THE 20-BAY TRUSS 


Weight, 


Design variables, 


Active constraints 

CPU 

Sensitivity 

lbf 



in. 2 





time, 



A-i 

Ai 

^3 

A 4 

^5 

Stress 

Displacement 

sec 


2023.2 

3.44 

7.03 

0.41 

1.38 

1.38 

3 

1 

2.5 

Determinate 

2021.8 

3.44 

6.99 

0.41 

1.39 

1.39 

3 

1 

14.8 

Analytical 


The SQP algorithm converged to the same solution for the determinate as well as for the analytical 
sensitivities. With the approximation, the CPU time to solution was 590 percent faster. With five design 
variables and four active constraints, this was a well-posed mathematical programming problem. 

Computational Efficiency 

All six examples converged to the correct solution with the approximate sensitivity. To examine 
the computational efficiency when the approximate sensitivity was used, we solved the three larger 
problems: the ring, the wing, and the 20-bay truss in a controlled environment using the SQP algorithm 
on an SGI workstation with an IRIX 6.5 operating system. The CPU time to optimum solution was 
measured for both the determinate and the closed-form analytical sensitivities. The CPU time to solution 
is depicted in table VII. 


TABLE VII.— CPU TIME TO SOLUTION FOR LARGE PROBLEMS WITH THE SQP ALGORITHM 


60-bar ring 

Forward-swept wing 

20-bay trass 

Sensitivity 

2.3 

8.7 

2.5 

Determinate 

7.4 

67.5 

14.8 

Analytical 


For the 60-bar trussed ring, the optimum solution was reached in 2.3 CPU sec with approximate 
sensitivity. The time increased threefold for analytical sensitivity. For the forward-swept wing, the 
time factor in favor of the approximation was almost eightfold; 8.7 CPU sec with approximation against 
67.5 sec for analytical sensitivity. The time ratio was 6 for the 20-bay truss; 2.5 and 14.8 CPU sec with 
approximate and analytical sensitivities, respectively. Overall, the time to solution increased from 
threefold to eightfold for the analytical sensitivity. Approximate sensitivity increased computational 
efficiency by several orders of magnitude. 


Convergence Pattern 

The convergence of weight versus CPU time to solution with analytical and approximate sensitivities 
for the three large problems is depicted in figure 8. For each problem, optimization was begun with the 
same initial design. Both methods produced similar optimum solutions. Consider the graph in figure 8(a) 
for the 20-bay truss. The convergence patterns, with and without approximation, portray undulations that 
are quite similar. However, the convergence is very rapid with the determinate sensitivity. The conver- 
gence characteristic is similar for the 60-bar trussed ring and the forward-swept wing shown in figures 
8(b) and (c), respectively. 
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Figure 8. — CPU solution time for SQP algorithm, (a) Twenty-bay truss, 
(b) Sixty-bar trussed ring, (c) Forward-swept wing. 
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Angle Between Search Directions 

The iterative optimization process moved along a search direction (y7) from one design point to 
another. The search directions, generated from the gradients, differed for the analytical | d anl j and 
determinate sensitivities [ d 'determinate ) • If the difference between d an j and ^determinate was small, then the 
contribution to sensitivity from the indeterminate factor could be considered to be negligible. Otherwise, 
the contribution from the indeterminate factor could be significant. To examine this issue, we defined an 
angle 9, at the i iteration between search directions generated using determinate ^determinate an d analytical 
d an i sensitivities as 


9 ; - = cos 1 

VI II I /; 

The angle would be zero (9, = 9) if the determinate and the analytical gradients were identical; 
otherwise, it would be nonzero (9, ^ 9). In some scale, the angle 9, was a measure of the difference 
between the closed-form and determinate gradients of the active constraints. For the five-bar truss shown 
in figure 3, the angle 9, was calculated at the initial and the optimum design points for the SQP algorithm. 
Numerical values for the directions and angle are given in table VIII. 


d an l ' d determinate 
^anl I d determinate I 


(19) 


TABLE VIII.— ANGLE BETWEEN SEARCH DIRECTIONS SOLVED 


WITH SQP ALGORITHM FOR THE FIVE-BAR TRUSS 


Direction vector, d 


Initial design point, 




Optimum solution point, 




angle Owtiai = 25° 




angle 0 ODt = 64° 



\ Sensitivity 

Sensitivity I 

I Determinate 

Analytical 

Detenninate 

Analytical | 


' 0.38' 



' 0.311 



1.51' 



-1.49 



-0.16 



-0.03 



-1.44 



0.19 


initial _ J 

^determinate — 1 

-0.23 

> 

j initial 
“anl - * 

-0.36 


i 0pt - J 

“determinate 

4.14 

>10“ 5 

( anl-opt 

{4} 

1.77 

•10“ 5 


-0.39 



-0.22 



-0.30 



-0.30 



-0.09 



-O.OVj 



0.00 



0.00 



The angle was different at the initial as well as at the optimal design points with values 9 in i ti ai = 25° 
and 9 opt = 64°. For the three large problems, the angles were calculated for the SQP algorithm and are 
given in table IX. At the initial point, the angles varied from 29° for the forward-swept wing to 5 1° for the 
29-bay truss. At the optimum solution point, the minimum and maximum values for angles were 39° and 
89°, respectively. In summary, the directions taken from the initial point to reach the optimum solu-tion 
were different for the analytical and the determinate sensitivities. In other words, the indeterminate 
component of the sensitivities in equation (4) was not negligible and it changed the path of optimization. 


TABLE IX.— ANGLE BETWEEN SEARCH DIRECTIONS FOR 
LARGE PROBLEMS SOLVED WITH SQP ALGORITHM 


Problem 

Angle, 0, deg | 

Initial design point, 

® initial 

Optimum 
solution point, 

0oot 

60-bar trussed ring 

47 

39 

Forward-swept wing 

20 

89 

20-bay truss 

51 

40 
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Discussion 


Design optimization was not sensitive to the precision of the gradients of the stress and displacement 
constraints of an indeterminate truss. Simple determinate sensitivities performed very well for all design 
problems. The CPU time to optimum solution was substantially reduced with the determinate sensitivity. 
The question is: Why did the determinate sensitivity outperform the full closed-form gradient? An answer 
is attempted through a discussion of the nature of the structural design problem. 

Nature of the Structural Design Problem 

The nature of the design optimization problem was examined by considering the three-bar truss 
shown in figure 2. The problem had three design variables (A i, A 2 , A 3 ). It had three stresses (gi, g 2, G3) 
and three stress constraints (g al , g a2 , g a i). Likewise, there were two displacements (X u X 2 ) and two 
constraints (g xl , g x2 ) for each load case. There were three implicit structural analysis relationships between 
the five behavior variables: 


Gi - G 2 + G 3 = 0 

(Ha) 

X 1 = 4( a2 _2(71 ) 
E 

(Hb) 

x 2 = — (<T>) 

(He) 


Equation (1 la) is the compatibility condition (CC), and it is expressed in stresses. Equations (lib) 
and (lie) are the deformation displacement relations, also expressed in stress variables. It is important to 
observe that the three implicit relationships do not explicitly contain the design variable. In other words, 
the gradient of one stress can be expressed in terms of the gradient of other stresses. Such as for example, 


(Vg 2 = V (gi + G 3 )) and 




If the second stress g a2 and second displacement 


g x 2 constraint become active, then the coefficient matrix [Q] of the direction vector { d | in equation (1) 
will become singular and [Vg a2 ; Vg v2 ] will become rank deficit. Four possible singularity situations are 
listed in table X. In case 1, a member stress (g 2 ) is dependent on a single displacement (x 2 ). The rank of 
the sensitivity matrix was 1, but it was restored to 2 with the approximation. Consider case 3 with three 
active constraints. The approximation had a full ra nk of 3, whereas the analytical sensitivity had a rank 
deficit at 2. The other two cases also exhibited deficient rank. 


TABLE X.— RANK OF SENSITIVITY MATRIX FOR A THREE-BAR TRUSS 


Case 

Set of active constraints 

Rank of sensitivity matrix | 


Total 

Number of each 

Constraints 

Analytical 

Detenninate 

1 

2 

1 stress and 
1 displacement 

So2i £x2 

1 

2 

2 

3 

3 stresses 

So 1? So2i So3 

2 

3 

3 

3 

2 stresses and 
1 displacement 

Soli So2i Sxl 

2 

3 

4 

5 

3 stresses and 
2 displacements 

Soli So2i So3i Sx2i Sx3 

2 

3 
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The four singularity cases listed in table X cannot be ascertained prior to the initiation of the 
optimization calculations because constraint activity depends on the value of design variables. Singularity 
will be avoided when determinate sensitivity is used because it restores the full ra nk of 
the sensitivity matrix as shown in table X. In other words, the implicit relationship of the behavior 
constraints induced singularity in the design optimization of this truss. 

Coefficients in Equilibrium and Compatibility Matrices 

The number of stress and displacement components in the implicit relationship, similar to that in 
equation (11), depends on two quantities: q ee and q cc . The number of entries (or nonzero coefficients) in a 
column of the equilibrium matrix [5] is q ee . Likewise, the number of nonzero coefficients in a row of the 
compatibility matrix [C] is q cc . Both q ee and q cc are small numbers. In other words, a small number of 
stresses are dependent. Likewise, a displacement is dependent on few stresses. Consider the example of 
the 20-bay truss shown in figure 7. The number of coefficients in a column of its equilibrium matrix 
varies between one and four (1 <q ee < 4; four is more prevalent). The coefficients in a row of the 
compatibility matrix range between 6 and 20 (6 < q cc < 20; the typical number is six). Six typical 
dependence relationships for the truss are given in table XL 


TABLE XI.— RANK OF SENSITIVITY MATRIX FOR A 20-BAY TRUSS 


Case 

Set of constraints 

Rank of sensitivity matrix 


Number 

Constraints 

Determinate 

Analytical 

1 

2 

<?ob 8x2 

1 

2 

2 

3 

8(54, 8x5, 8 x6 

2 

3 

3 

5 

8(5 68? 8x69,8x7 0 ? 8x7 U 8x72 

4 

5 

4 

6 

8g l, 8 g2, •••? 8o 6 

5 

6 

5 

6 

8o3l, 8(5 32, •••> 8(5 36 

5 

6 

6 

20 

8(5 2, 8o7, 806 7 

19 

20 


For case 1, one stress and one displacement are dependent because q ee = 1. For case 2, one stress 
is dependent on two displacements because q ee = 2. For case 3, one stress is dependent on four dis- 
placements because q ee = 4. For cases 4 and 5, six stresses are dependent because q cc = 6. For case 6, 

20 stresses are dependent because q cc = 20. The sixth case is interesting because the 20 stresses belong to 
the bottom chord members of the truss. The bottom chord is a natural load path but it can promote singu- 
larity in the optimization process, which however, can be avoided when approximate sensitivity is used. 
For a truss, a small number (in the range of two to six) of active stress and displacement constraints can 
be dependent. 

A traditional structural optimization problem contains dependent constraints. A small number of 
active stress and displacement constraints can be dependent. The multitude of implicit constraints reduces 
the rank of the coefficient matrix \Q\. Simple determinate sensitivity worked well because it restored the 
full ra nk for each of the six cases shown in table XI. Earlier, we suggested (ref. 9) that a set of 
independent constraints should be separated out of the given stress and displacement constraints by using 
a singular value decomposition algorithm. The exercise has to be performed before the generation of each 
search direction. This technique works well for small problems. For larger problems, the decom-position 
process increases the numerical burden in optimization, which is already computationally intensive. The 
current recommendation is to use simple determinate sensitivity because it restores the full rank of the 
matrix [ Q\ and converges more rapidly. 
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Adjustment for the Stiffness Method 

Design is stress driven both for stress and displacement limitations. The area of a truss bar for stress 

limitation can be updated as A new = — .4° 1(1 . Likewise, the displacement formula ( X = JGF) can be 

<*0 

manipulated to obtain an area update formula for the stiffness limitation (ref. 13). The two features make 
the method of force an attractive tool for design applications. We, however, realize that the stiffness 
method is very popular. The question is, “Can sensitivities be approximated when the stiffness method 
is used as the analysis tool in optimization?” Such an approximation is straightforward for the stress 
constraints. It may pose a challenge for the displacement limitation because it is a global variable. In the 
stiffness method, the stress sensitivity can be obtained by dividing the force or stress parameters (F or a) 
by the square of area or area, respectively, -FI A 2 or -a! A. The force or stress output of a stiffness code 
can be adjusted to obtain the approximate sensitivity for the stress constraints. 

The difficulty encountered in calculating the approximate displacement sensitivity is illustrated 
by considering the three-bar truss as an example. Consider one term in the derivative of the first 
displacement X\ with respect to the area A\. For simplicity, one load component is set to zero, P v = 0. The 
closed-form derivative can be written as 


6X, 

dA 1 


2 (Vi A\ + 4A 2 A 3 + 2 42 Al \ IP X 

\ 2 

E (42A 2 A 3 + 2A\A 3 + 42A x A 2 j 


( 12 ) 


The procedure of separating the sensitivity into determinate and indeterminate factors in equation (8) 
cannot be directly extended to the stiffness method, such as for example in equation (12). The stiffness 
method, in general, appears to have two major impediments for design calculations. 


(1) The method has fewer equations m than the number of design variables n: m < n. The three-bar 
truss has three design variables, three bar stresses, but there are only two stiffness equations. Three 


equations are required to size the three bar areas, like A\ = 


F}_ 


, A 9 — 


Vi 

<*o 


and A 3 - 


3 . 


. In other words. 


it is not easy to link bar areas to displacements. At best, this link would provide a relationship of three 
design variables to two displacements, and it would not be a one-to-one mapping. 

(2) The stiffness method does not allow free movement between analysis variables like 1FM, which 
allows movement from force to displacement, { X] = [.7][G] { F}, and vice versa, {F} = [G] '[B] r {X\. The 
two formulas, along with their governing equation, [5] {F} = \P \ , make 1FM a very attractive tool for 
sensitivity calculation and design optimization. 


In the stiffness method, the sensitivity expression given by equation (9) can be used provided the 

matrix [.7] can be approximated. For static response only, it can be approximated as [j] = 

The generation of the equilibrium matrix is straightforward because this, in essence, is the concatenation 
of the transformation submatrices used to change the local to the global coordinate systems. The inverse 
of [ BB t ] cannot be avoided, except that [B~\ is a very sparse matrix. 

Extension to Other Structure Types 

Extension of the approximate expressions is straightforward for beams and framework. Consider, for 
example, a beam with moment as the analysis variable and moment of inertia / and depth d as the design 


[BBT ]- 1 
„ [ 0 ] 
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variables. The flexure formula can be differentiated to obtain the approximate sensitivity for stress 
constraints: 


G = 



da 

~dl 


Md 
2 1 1 


■ + 


neglect — »■ 


d 8M^ 
21 dl , 


( 13 ) 


The derivative of the moment in the displacement formula (X=JGF) would provide the approximate 
sensitivity of the stiffness constraints. In other words, in 1FM, the generation of the closed-form and 
approximate expressions for beams is quite straightforward. The logic can be extended for framework that 
would require a mixing of the truss and beam expressions. Then computer software has to be developed to 
compare design optimizations using approximate and analytical sensitivities for such structures. The 
exercise is worth the effort because singularity can be eliminated to make optimization robust for flexural 
structures. 


Conclusions 

There are numerous dependent constraints in the design optimization problem of an indeterminate 
truss. A small set of stresses can be dependent. A stress also can depend on a few displacements. A truss 
design with many sets of dependent constraints may not be a well-posed mathematical programming 
optimization problem. However, it is a real-life industrial design problem. In optimization calculations, all 
constraints should be used in defining the feasible region. Sensitivities of only the independent constraints 
should be used to calculate the direction vector. The independence criterion will be satisfied when the 
proposed simple determinate design sensitivity is used. The optimum solution was reached with the 
determinate sensitivity even though the search directions differed for the determinate and analytical 
sensitivities. The use of simple determinate sensitivity substantially reduced the CPU time to solution. 

The integrated force method is an efficient analysis tool for the calculation of determinate sensitivity in 
particular and for design application in general. The concept of using approximate sensitivity in design 
optimization should be extended to flexural structures like beams and framework. 
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Appendix — Symbols 


A i bar areas 

B ] m x n equilibrium matrix 

[C] r x 77 compatibility matrix 

d search direction 


^determinate > ^anl 

E 

F 

[G] 

gt 

[Vg] 

[J] 

m 

n 

{ P } 

Qcc 

Qee 

r = n — m 

[s\ 

Uj 

Ujo 

w 

x, y,z 


Pi 


Oi 


Oo 


search directions for determinate and analytical sensitivities 

Young’s modulus 

member force 

77 x 77 flexibility matrix 

7 th constraint 

sensitivity matrix 

first 777 rows of [[5] '] r 

number of displacements 

number of internal forces 

load vector 

number of entries in a row of the compatibility matrix [C] 

number of entries in a column of the equilibrium matrix [5] 

number of dependent members or compatibility conditions 

77 x 77 IFM governing matrix 

nodal displacement 

limitation on nodal displacement 

weight 

coordinates 

weight density for 7 th bar 
bar stress for 7 th bar 
allowable stress 
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